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Abstract 



The pionic decay of a possible d' -dibaryon in the process d' — > N + N + tt is 
studied in the microscopic quark shell model and with a single-quark transition operator 
describing the transition q — > q + tt. For the d' with quantum numbers J p = CP, T = 0, 
we employ a six-quark shell-model wave function with a spatial s 5 p[51]x-configuration 
' and with N = 1 harmonic oscillator quanta. It is shown that the pionic decay width 

\& . depends strongly on the mass and size of the d' . In the case that the calculated d' 

mass is close to the experimental one a small pionic decay width of l\ w 0.04 MeV 
is obtained. This is an order of magnitude smaller than the experimental T° xp ~ 0.5 
MeV. Two possibilities to improve the calculated width are suggested. The effect of the 
nonstatic correction term in the transition operator and the influence of the form factor 
' at the decay vertex on the decay width are also discussed. 



1. INTRODUCTION 

Dibaryons have a long history of more than three decades. The experimental search for 
dibaryons in the nucleon-nucleon (NN) channel has however not been successful yet. On 
the other hand, more exotic baryon-number B = 2 systems, that is dibaryons which have a 
genuine six-quark (n(q) — n(q) = 6) structure that cannot be described by hadronic degrees 
of freedom are predicted by quantum chromodynamics (QCD). If such an exotic dibaryon is 
proven to exist in even a single case, the consequences would be far reaching for understanding 
QCD in the low-energy domain of hadron and nuclear physics Q. 

Recently, new pionic double charge exchange (DCX) (7r + ,7r~) data have revived the in- 
terest in dibaryon physics. DCX measurements on a number of nuclear targets ranging from 
4 He to 56 Fe show a sharp peak in the excitation function at an incident pion energy T n ~ 50 
MeV and at forward angles #l ~ 5° [||, [|. The location of the peak at T w ~ 50 MeV is 
nearly independent of the nuclear target. In the 4 He - target case, the peak appears around 
T n ~ 80 — 90 MeV Q due to the additional 28 MeV of binding energy required to break up 
the 4 He nucleus into an unbound four-neutron final state. 
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The DCX process necessarily involves two nucleons due to charge conservation and is 
sensitive to short-range NN correlations. In addition, because the peak position and the 
width are largely independent of the target, the observed resonance-like peak is likely to be 
connected with some elementary process. The invariant mass spectrum of the p + p + tt~ 
subsystem in the p + p^p + p + ir~ J t-ir + reaction also shows a peak at the same position 
as the (7r + , 7T _ ) DCX experiments, which supports the interpretation of the observed peak as 
an elementary process 

The resonance peak, which has been called d' -dibaryon, has an energy = 2065 
MeV, spin-parity J p = 0~ and isospin T = ||]. The width of the peak is reported to 
be T me dium ~ 5 MeV when the effect of the Fermi motion is subtracted Q. Such quantum 
numbers and such a small width are incompatible with the picture of an ordinary baryon- 
baryon ( NN*, N*N*, etc. ) bound state. 

Theoretically Mulders et al. predicted that six-quark states with quantum numbers J p = 
0~,T = should have a tetraquark (q 4 ) — diquark (q 2 ) stretched structure and a mass 
M ~ 2100 MeV within the MIT bag model ||. Recently, the Tubingen group has studied the 
microscopic structure of a possible d' -dibaryon in the nonrelativistic quark model (NRQM). 
Two different bases have been employed, the quark cluster model and the quark shell model 
|], 0] in both cases paying due attention to the Pauli principle. In ref. || it has been shown 
that the nonrelativistic quark cluster model with a q 4 — q 2 -cluster configuration can account 
for the d' -mass of M# ~ 2100 MeV if the size parameter of the d' wave function is properly 
chosen ||. In the translationally invariant quark shell model (TISM) |], 0] the six-quark 
states with J p = 0~,T = are evaluated in a model space including N = 1 and N = 3 
harmonic oscillator (h.o.) excitations. The same effective interaction potentials as in the 
cluster model approach are used. The quark shell model gives similar d' -masses as the 
quark cluster model. 

Another interesting problem is whether the NRQM can explain the observed decay width 
of the d' -resonance. This is the objective of this article. It is reported that the major part 
of the observed width T me( num in DCX is due to the spreading width. The spreading width 
is a result of the interaction of the d' with the nucleons in the nucleus; it is absent if the 
d' is formed in an elementary reaction such as p + p — > p + p + vr~ + ir + . Since the d' is 
not allowed to decay into the NN channel because of its quantum numbers J p = 0~, T = 0, 
the dominant decay mode of the d' would be the NN7T channel. The pionic decay width of 
the d' is estimated as small as T n ~ 0.5 MeV [||. We evaluate the pionic decay width by 
employing the microscopic quark shell-model wave function of the d' as obtained in Ref. J7j. 

In contrast to the isospin-0 assumption for the d' used here, Valcarce et al. || have 
discussed the possibility that the J p = 0~ resonance in the NN-zr system seen in DCX 
experiments is an isospin-2 resonance in the nucleon-A channel. Using a nonrelativistic quark 
model they have shown that the nucleon-A interaction is attractive in the J p = 0~ channel 
and that the decay width is explained if the calculated mass of the resonance approaches the 
experimental one. 

This paper is organized as follows. In sect. 2 the basic expression for the pionic decay 
width is derived and its approximations are discussed. A refined treatment of the transition 
operator is also given in sect. 2. The quark shell model wave functions and the calculated 
d' -dibaryon mass are given in sect. 3. In sect. 4 the calculated results of the pionic decay 
widths are shown and the relation between the width and the d' -mass are discussed for 
several parameter sets. A summary is given in sect. 5. 
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2. PIONIC DECAY WIDTH OF THE d' -DIBARYON 



2. 1 Basic formula 

The dominant decay mode of the d' -dibaryon with quantum numbers J p = CP , T = 
is the pionic decay, such as d' — > p + p + ir~ , p + n + 7r° and n + n + ir + . We treat the 
decay operator as well as the dibaryon and final nucleon states in terms of quark degrees of 
freedom. The decay of the d' -dibaryon arises from the single-quark transition q — > q + n, 
which is depicted in Fig. 1. 



Fig. 1 



Fig. 1 The single-quark q — > q + tt transition process, p (p') denotes the initial (final) quark 
momentum and k is the pion momentum in the final state. 



The transition operator is expressed after nonrelativistic reduction as 
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° = Ztj«rmTy^--^-^ 7 ^. (1) 

Here, f nq represents the coupling constant at the qq7r- vertex, E n (m n ) the pion energy (pion 
mass) and c() the isovector of the pion field. The <t(t) represent the spin (isospin) operator 
of a single quark. 

The nucleons and the pion emitted in the d' decay d' — > N + N + ir have momenta 
q 1} q 2 and k n , respectively, in the CM. frame. We define 

Q = qi + q 2 , 9=^"^ (2) 

where Q(q) is the total (relative) momentum of the two nucleons. Momentum and energy 
conservation of the decay process are expressed as 

= Qi + q 2 + K = Q + K , (3) 

respectively, where represents the d' -dibaryon mass. 
The pionic decay width is expressed as 

rr = 2. / dQdqdKS (M d , - 2M N - - - VW + kl) 

1 Mf Mi Tf 

x | < " NN" JfMfTfTf ; J l M l T l r l > | 2 

= 2vr / dqdK5 U, -^ N -^-^-^ml + kA 
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x o 7 + 1 2 S XI 

x | < " NN" JfMfTfTf, w |O|* d ', > | 2 



(5) 



where the initial d' quantum numbers are referred to as J{M{TiTi and "NN" refers to 
the final two-nucleon state with quantum numbers JfMfTfTf. To be specific we assume 
j/ 3 = CT,Tj = and jj 3 = + ,Tj = 1 due to the isovector nature (r) of the transition 
operator O. Then the total pionic decay width can be written as 

V = IV + r> + rv (6) 

and IV is expressed as 

IV = 2vr / dgdM (M d , 21/, J- \ • fc 2 ) 

x | < "NN"J f = + ,T f = l,r f = 1;<V \0 \ 

* df ,J i = 0-,T i =0,T i = 0>\ 2 . (7) 

Similar expressions hold for IV and T w + . Owing to isospin invariance, one has 

IV = r 7T° = V + (8) 



5.2 Restriction to dominant intermediate six- quark states 

The decay transition matrix element can be evaluated by inserting a complete set of 
six-quark states as 

< "AW" J f = + ,T f = l,Tf = l;<l> v -\O\* <e ,Ji = 0-,T i = 0,Ti = 0> 
= ^<"iViV"J / = + ,T / = l,r / = l|V6 g ,n> 

n 

x<il> 6qtn ,<f> 1[ -\O\V d ,,J i = 0-,T i = 0,T i = 0> . (9) 

The first factor on the r.h.s. of Eq. (9) is the overlap between the antisymmetrized 
two-nucleon state and the intermediate six-quark state of the n-th excitation. The second 
factor which contains the decay dynamics connects the d' -dibaryon and the intermediate 
six-quark state accompanying the emitted 7r~-meson wave. The quark shell model s 5 p[51]x- 
configuration with N = 1 h.o. excitation is used for \Pd'(J = 0~, T = 0). 

The following arguments restrict the sum over intermediate states in Eq. (9). First, the 
pionic decay proceeds through a one-body operator. This means that only such six-quark 
configurations in which a single orbit (s or p) is changed from the initial s 5 p-configuration 
are allowed. Second, those intermediate six-quark states which have a large overlap with 
asymptotic AW components (large fractional parentage factors) will give an appreciable 
contribution to the transition matrix element. In view of this, the two lowest-lying six-quark 
states are adopted as intermediate states in this work. These are the Ohu 4>Q q {s Q [Q]x ,J = 
0+, T = 1) six-quark state with N = and the 2 hu) excited state ipQ q (s 4 p 2 [42] x , J = + ,T = 
1) with N = 2. A more detailed discussion of these wave functions as well as of the d' wave 
function is given in sect. 3. 
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2.3 Nonstatic correction and form factor in the transition operator 



One can go beyond the static limit of Eq. (1) and consider the next-to-leading order 
relativistic correction in the transition operator O. The nonstatic correction term at the 
qq7r-vertex in Fig. 1 might be important in the present quark model description of the d' 
decay because the ratio of pion and quark masses is not small. The nonstatic correction of 
order 0( 1 ^ L ) can be taken into account by modifying the momentum-dependent part of the 
operator O of Eq. (1) as follows 



(<Tj ■ k) 



(k 



m 7 



P' +Pn 



2m n 



(10) 



where p (p 1 ) is the initial (final) quark momentum ||. The quark momentum-dependent 
term in Eq. (10) is often referred to as the recoil correction. 

Another refinement is due the finite size of the qq7r- vertex in Fig. 1. Usually, this is 
treated by defining a renormalized momentum-dependent coupling in place of the f nq in Eq. 

(1) 



fnq y f- 



A 2 



.1/2 



A 2 + fc 2 



(11) 



consistently with the pion-exchange potential [1C]. The cut-off mass A measures the extension 
of the core-size of the quark-pion interaction. 



3. THE d' -DIBARYON MASS AND QUARK SHELL MODEL WAVE FUNCTIONS 

3.1 The d -dibaryon mass and wave function 

The d' -dibaryon mass and wave functions which have been studied in Refs. H are 
employed in the present calculation of the pionic decay rate of the d' . The basic ingredients 
and the results of Refs. || 0] which are essential for the present calculation are thus briefly 
recapitulated. 

Although quark dynamics is known to be described by QCD, the many-quark system 
(n(q) > 3) has not been solved successfully from first principles except for calculations in the 
lattice gauge theory. In the low-energy domain, the nonrelativistic constituent quark model 
works well in accounting for the regularities of single hadron spectra. 

The translationally invariant quark shell model Hamiltonian which includes the effective 
quark-quark interaction is written as 

n T) 2 P 2 n n 

H = EK + ^-) - + E v Cmf (.^ T,) + E v^iruTj). (12) 

The confinement potential V Con ^ is necessary; it is assumed to be linear @ or quadratic H 
in |r| = | ri — TVj|, i.e., 

V Con f{r urj ) = -l^.X^aJn-rjl-C), (13) 
or 

= -a' c X i -\ j {r l -r j ) 2 , (14) 
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where Aj is the SU(3) co ; or matrix of the i-th quark. The residual interaction V res is composed 
of the chiral field ( it and er ) potential and the one-gluon exchange potential. The former 
one-pion and one-sigma exchange potentials are closely related to the spontaneous breaking 
of chiral symmetry of QCD [|ll], |1^| . The latter one-gluon exchange potential is responsible 
for the short-range part of the interaction [jl^]. They are expressed as 



^ (ri , n) = |^_^ (Tj . Tj)(<Ti . Vr)(<Tj .. Vr) (£^_f^) , (15) 



with 



4lT 47T 



(2m q ) 2 + ml, A n = A a = A, 



(17) 



and 



V OGEP {r urj ) = -asX, ■ A J - - -^(1 + \a t ■ a^r) \ , (18) 
4 I r mf: 6 I 

where <Xj is a spin operator. The vertex form factor ( A2 A ^ ) 1 ^ 2 is used in deriving Eqs. 

(15) and (16). The relations g 2 q = f 2 q {^^-) 2 and f nq = ^f w N hold, where f n N is the NN7T 
coupling constant. The parameter sets used in this work are listed in Table I for both linear 
and quadratic confinement potentials. For sets 1-2 the parameters are determined by fitting 
the light baryon masses ( N, A, N*(1535) ), and by satisfying the nucleon stability condition 
dM]\r I 'dbfq = and the approximate constraints for A and m q . Parameter sets 3-5 fit the 
N and A masses but underestimate the N*(1535) mass by about 300 MeV. The latter 3 
parameter sets all have a small confinement strength a c . 



Table I Quark shell model parameters and the calculated d' -dibaryon mass M^- The 
6jv and are the h.o. size parameters of the final nucleon and the d' -dibaryon wave 
functions, respectively. m q is the constituent quark mass. a s ,a c , C, a' c and A are the potential 
parameters. Set 1-3: The values of the parameters and the corresponding d' -dibaryon masses 
for a linear confinement potential (Eq. (13)). Set 4-5: The values of the parameters and the 
corresponding d' -dibaryon masses for a quadratic confinement potential (Eq. (14)). 





b N 


a s 


a c 


C 


m q 


A 




M d , 


SET 


[fm] 




[MeV • fm" 1 ] 


[MeV] 


[MeV] 


[fm" 1 ] 


[fm] 


[MeV] 


1 


0.45 


0.127 


462 


549 


338 


5.07 


0.59 


2705 


2 


0.47 


0.074 


423 


500 


296 


7.60 


0.65 


2680 


3 


0.60 


0.816 


25 


-32 


313 


10.14 


1.24 


2162 




b N 


a s 


a' c 




m q 


A 


b 6 


M d , 


SET 


[fm] 




[MeV • fm" 2 ] 




[MeV] 


[fm" 1 ] 


[fm] 


[MeV] 


4 


0.595 


0.958 


13.97 




313 


4.20 


0.78 


2484 


5 


0.595 


0.958 


5.0 




313 


4.20 


0.95 


2112 
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The lowest-lying harmonic oscillator state which is compatible with d' quantum numbers 
J p = 0~ and T = involves one quantum excitation N = 1 and is uniquely expressed as 

y d ,(N = l,s 5 p[5l] x ,(\n) = (W), L = l, S = 1,J = 0~,T = 0, 

[2211] CT , [222] c ). (19) 

Here, [...]x> [•••]ct and [...]c specify the space-, color-isospin- and color- symmetry, respec- 
tively, and (XfJ,) represents the Elliott SU(3) label of h.o. quanta. The h.o. size parameter of 
the *$>d' wave function is b^. Other higher excited configurations with N = 3 excitations and 
[42] x symmetry are considered in Ref. |/]]. A lower d' -mass is obtained in the configuration 
mixing calculation but the mixing amplitudes are small. Therefore, in the present paper, we 
adopt only the lowest configuration of Eq. (19) for the d' -wave function. The size parameter 
&6 is determined by minimizing the d' -mass. Calculated d' -masses and corresponding 
size parameters b§ are shown in Table I. It is noticed that in the cases of weaker confinement 
potentials (set 3 &; 5), the size parameters b% become large and consequently lower dibaryon 
masses M$ , which are comparable with the experimental resonance energy of 2065 MeV, are 
obtained. 

3.2 The intermediate six- quark wave functions and the final two-nucleon state 

Among all possible intermediate states only the lowest two six-quark states with J p = 
+ , T = 1 are adopted in Eq. (9) for the evaluation of the pionic decay width of the d' . The 
criteria for selecting these six-quark states have already been mentioned in sect. 2.2. The 
lowest two states with N = and N = 2 h.o. excitations are expressed as 

V^i = <M N = 0,s 6 [6] x ,M = (00), L = 0, 5 = 0, J = 0+,T = 1, 

[222] CT , [222] c ) (20) 

and 



^6,,2 = iM A r = 2,sV[42]x,(A / u) = (20), L = 0, S = 0, J = 0+T = 1, 

[42] ct, [222] c ) ■ 



(21) 



The probability amplitudes for finding the two-nucleon component in these two six-quark 
states (fractional parentage coefficients) are y/1/9 and — a/1/25, respectively [14, if a 
common h.o. size parameter is assumed for ipeq,i(i = 1>2) and if the nucleon has an s 3 - 
configuration. 

The final outgoing two-nucleon state is approximately written as 



| "iViV"; J =0+,T=l> 

L = 0,5 = 0, J = 0+,T = 1 > . (22) 

In Eq. (22), an exchange of three-quarks between nucleons is taken into account. ipN repre- 
sents a free nucleon and its size parameter is given by b^ in Table I. 
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4. RESULTS AND DISCUSSION 



The expression for the pionic decay width of the d' is listed in the Appendix. The analytic 
form of this expression is an advantage of the present model; it facilitates the discussion and 
interpretation of our results. The decay widths r„— and T n are evaluated by adopting the 
wave function vl/^' of Eq. (19) and several sets of parameters for the linear confinement (sets 
1 — 3) and the quadratic confinement (sets 4 & 5) potential cases. The calculated decay 
widths are shown in Table II under the abbreviated symbols S, S+NSC and S+NSC+FF, 
where S refers to a calculation employing the static decay operator of Eq. (1), NSC includes 
the nonstatic correction term of Eq. (10), and FF refers to a calculation using the form 
factor of Eq. (11). In Table II, two numbers for T n - are listed in each column S, S+NSC 
and S+NSC+FF in order to show the effect of the intermediate six-quark states in Eq. (9); 
the numbers without parenthesis are widths evaluated by adopting both configurations t^eg,! 
and ip6q,2 (abbreviated as I & II, respectively) and the numbers with parenthesis [ ] are the 
widths evaluated with configuration I only. 

Table II Calculated 7r _ -decay width r„.- and the total 7r-decay width T n of the d' -dibaryon 
for the adopted parameters sets 1-5. The experimental data (2| is also shown. The symbol 
S denotes the calculation with use of the operator O of Eq. (1) in the static approximation, 
S+NSC the calculation with the operator O with inclusion of the nonstatic correction term 
of Eq. (10) and S+NSC+FF the calculation with the operator O with inclusion of nonstatic 
correction term and the form factor of Eq. (11). The numbers without the parenthesis in each 
column of S, S+NSC and S+NSC+FF of r^- are the calculations when the configurations 
I and II (ipeq t i and ^69,2) are adopted for the intermediate six-quark states. The numbers 
with the parenthesis [ ] are the calculations when only the configuration I is adopted for 
the intermediate six-quark state. The right-most column lists the calculated decay width 



SET 


6 6 [fm] 


M d ,[MeV] 




IV [MeV] 




r w [MeV] 


S 


S+NSC 


S+NSC+FF 


S+NSC+FF 


1 


0.59 


2705 


8.16 [8.83] 


7.16 [7.85] 


5.33 [5.82] 


16.0 


2 


0.65 


2680 


6.70 [7.15] 


6.34 [6.84] 


5.53 [5.95] 


16.6 


3 


1.24 


2162 


0.03 [0.03] 


0.01 [0.02] 


0.01 [0.02] 


0.04 


4 


0.78 


2484 


4.22 [4.52] 


3.74 [4.06] 


2.90 [3.13] 


8.69 


5 


0.95 


2112 


0.02 [0.02] 


0.01 [0.01] 


0.01 [0.01] 


0.04 


EXP. 












~ 0.5 



Let us first discuss the decay width calculated with the operator O in the static approx- 
imation in Eq. (1). The calculated numbers are quoted in column S in Table II. The decay 
amplitude which connects the dibaryon state ^d' an d the final two-nucleon and one-pion 
state is the essential ingredient in determining the pion decay width. The comparison of the 
calculated decay width with experiment provides a severe test of the validity of various model 
assumptions. 

The pion decay width is very sensitive to the d' -dibaryon mass which essentially 
determines the phase space of the three-body NN-7r-decay. The decay kinematics is given in 
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Eqs. (3) and (4). Due to the larger phase space, a large M$ leads to a large pionic decay 
width r„.- and IV This is seen in the cases of sets 1 & 2 (linear confinement) and set 4 
(quadratic confinement) in Table II. If the calculated M^i is close to the experimental value 
2065 MeV as in set 3 (linear confinement) and set 5 (quadratic confinement), a small IV- 
(IV) is obtained. 

Another factor which affects the decay width is the overlap factor between the final two- 
nucleon state in free space and the two three-quark-states (s 3 ) in the six-quark configuration. 
The overlap-squared is expressed as 

.i + (WMV ' [ } 

where be(b]sr) is the h.o. size parameter of the six-quark (free-nucleon) state. The overlap- 
squared values for sets 1, 2 and 4 are almost similar and as large as ~ 0.6, while for sets 3 & 
5 the values are as small as 0.05 and 0.28, respectively. The latter small values for sets 3 & 
5 are due to the very different sizes of b§ from b^ as seen in Table I. 

In view of the above two key-factors, large decay widths IV- are obtained for sets 1, 2 
and 4, while much smaller rv are obtained for sets 3 & 5 (see the column S in Table II). 

One can see the role of adding the configuration II from the comparison of the two numbers 
without and with parenthesis [ ]. The inclusion of configuration II reduces the decay width 
r„.- by about 7 — 10 % compared with the values for configuration I only. This reduction 
can be physically understood. It is due to the destructive interference between the s 6 and 
the s^p 2 six-quark component of the NN wave function at short distances which results in a 



smaller intermediate six-quark amplitude [16]. The effect of the excited N = 2 intermediate 
six-quark state ^69,2 with an s 4 p 2 [42] x -configuration turns out to be significant but does not 
change the results qualitatively. This is because the transition amplitude which connects 
with ^6g,2 and (j> n is small compared with that between ^>d' and ^69,1 and 4> n . 

We proceed to discuss the effect of the nonstatic correction (NSC) term of the transition 
operator to the width. Since in the preceeding paragragh the contribution of ^6g,2 to the 
amplitude is found to be rather small, we evaluate the NSC effect only for the transition from 

to V>6g,i and <j) n . 

The right-hand-side operator in Eq. (10) is rewritten as 

K' • k ) - 

"*w> ■-■&>■*>}■ (24) 

The first term (o - ,- • k) is the operator in the static approximation and the second term the 
nonstatic correction. The transition amplitudes between and ^69,1 and <\> v are expressed 
symbolically in the following form : For the operator (<r • k) it reads 

< ip 6g ,l, 4hr I («r • fe) ■ • • I ¥ d / >~ const • (b 6 k 2 ) e -^ {b(ik)2 , (25) 
while for the operator due to the nonstatic correction it reads 

I { ■*) -(^X* -P) 

„ const- {^(p + p'^kf^e-^r (20) 
2m q 6 6 
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where b§ is the size parameter of \lv and V^i- 

Apart from the exponential factor, the following features are noticed in Eqs. (25) and 
(26). In the domain of large k, the term proportional to k 2 in Eq. (25) contributes most, 
while in the small k region the nonstatic correction dominates over the static approximation 
term due to the presence of the fc-independent term in Eq. (26). 

Combining Eqs. (25) with (26), the effect of incorporating the nonstatic correction term 
in the operator is written in a factorized form as [NSCI.]. An explicit expression for [NSCL] 
is given in the Appendix. Note that k is replaced by fco in the Appendix. The factor [NSCL] 
equals to 1 if no NSC is considered. The factor [NSCL] changes appreciably as k varies. The 
NSC reduces the decay width r„.- roughly by ~ 10 % for sets 1, 2 and 4 compared with the 
corresponding values in column S. However, for set 3 the NSC is quite important and the 
width r^- is reduced to half the value in column S, while for set 5 the NSC changes T w - only 
slightly. 

Finally, the form factor (FF) effect is discussed. The FF at the qq7r-vertex is of square- 
root type as shown in the Appendix. It is most effective for large momentum transfers k$. 
Thus a strong FF induced suppression of the decay width is seen for sets 1, 2 and 4 in which 
large momentum transfers fco are involved. On the contrary, the FF effect is negligible for 
sets 3 & 5 where the momentum transfer to the pion is small. 

The final results for the total pionic decay width are listed in the right-most column 
under the symbol S+NSC+FF, which are to be compared with the experimentally reported 
width T| xp ~ 0.5 MeV. The calculated T n for sets 1, 2 and 4 are too large in comparison with 
experiment. With and the parameters for sets 1, 2 and 4, the calculated dibaryon masses 
M^i are 600 — 400 MeV above the resonance mass 2065 MeV. This overestimation of the d' 
mass is mostly responsible for the large pionic decay widths. The calculated rV for sets 3 & 5 
are both 0.04 MeV which is an order of magnitude too small compared with the experimental 
data. Although the dibaryon masses close to the experimental one are obtained by using 
^d' an d the parameter sets 3 and/or 5, the corresponding sizes of the dibaryon b§ for these 
cases are large compared with the size parameter of the nucleon 6jv- Thus the overlap factor 
involved in the decay width calculation comes out too small for sets 3 &; 5 ( 0.05 &: 0.28, 
respectively ), which in turn results in a too small decay width. 

We now come to a critical discussion of these results. In order to explain the pionic decay 
width of the d' -dibaryon, first of all the experimental d' -mass has to be satisfactorily 
reproduced by the wave function vf^' . With the parameters of sets 1 , 2 and 4 which produce 
large d' -masses far from the experimental value it seems impossible to explain the decay 
width r^ xp . The d' mass can be explained by ^d' with parameter sets 3 & 5. It is however 
noted that very weak confinement potentials are adopted for these cases and the low d' 
-masses close to the experimental resonance energy are obtained at the price of having 
an extended radial d' wave function ^d' with a large b$. Thus, although the d' -mass M# 
is correctly reproduced for sets 3 & 5, the size-parameter mismatch between the final free 
nucleons (6jv) an d the three-quark-states (s 3 ) in ^>d' (be) results in a very small overlap and 
consequently leads to a small pionic decay width. 

Suppose that a certain Hamiltonian would produce the d' mass of 2112 MeV (set 5) with 
a size parameter b$ = 0.65 fm (0.595 fm) in the wave function ^/ . In this case one would 
obtain T T as 0.36 MeV (0.40 MeV) which is close to the experimental r° xp ~ 0.5 MeV. 

Two possibilities are suggested to achieve a small decay width. One possibility is to study 
the effect of additional quark-quark interactions which are not considered in the present work 
and which might explain the d' mass correctly. As we have argued above, the corresponding 
wave function ^d' should not be much wider than the wave function of a single free baryon 
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in order that the overlap factor of Eq.(|23[) remains close to unity. The second possibility is 
concerned with an improved treatment of the final state. In the present work the relative 
motion of the final state nucleons and of the pion are simply assumed to be described by 
plane waves. In a more complete treatment, the final state should be calculated with the 
same Hamiltonian that is used to calculate the mass of the d' ; in other words an NN-7T 
scattering wave function should be used in the calculation of the decay amplitude. If such 
final state wave functions are incorporated in the pionic decay process, for example in the cases 
of sets 3 & 5, an enhanced decay width would be obtained. This enhancement is expected 
because the relative wave functions are distorted and contain high momentum components 
which will in turn lead to an enhanced decay width P7| . It is noted that Schepkin et al. 
[ IS | have estimated the enhancement factor of the pionic decay width (77 ~ 5 in their model) 
when they considered the final state interaction for the decaying nucleons. 



5. SUMMARY 

The pionic decay width of a possible d' -dibaryon found in (7r + ,7r~) DCX reactions 
on many nuclear targets is studied in the microscopic quark shell model. We employ a 
single-quark transition operator which expresses the conversion of a single-quark state into 
another single-quark one by emitting a pion. We adopt the shell model wave functions 
devised by the Tubingen group. They have investigated the d' -mass and its structure in the 
translationally invariant quark shell model which includes the quark-quark interactions such 
as the confinement potential, the chiral field (n- and cr-exchange) and the one-gluon-exchange 
potential. In this work, the orbital d' -wave function with J p = 0~, T = is given by a 
s^p [5 l]x -configuration with N = 1 harmonic oscillator excitation. The pionic decay operator 
in its simplest form is obtained by the non-relativistic reduction of the qq7r-vertex and by 
imposing the static approximation on the operator. A more refined transition operator which 
includes the non-static (recoil) correction term and also the form factor at the qq7r-vertex is 
discussed. 

The pionic decay transition matrix element is evaluated in our model as the sum of 
products of two amplitudes; the first amplitude describes the transition from the d' -state 
to the intermediate six-quark state and a pion, the second amplitude represents the overlap 
between the intermediate six-quark state and the final two-nucleon state. The sum extends 
over important intermediate six-quark configurations. The calculated decay width turns out 
to be strongly dependent on the two basic properties of the d' , i.e. its mass M# and its 
size &6- The d' -dibaryon mass essentially determines the phase space of the pionic decay 
and hence controls the decay width. The size of d' strongly affects the overlap between the 
intermediate six-quark state and the final two-nucleon state. The pionic decay width also 
depends on the number of intermediate states included and on the treatment of the decay 
operator. 

The pionic decay widths and masses of the d' are evaluated for microscopic d' wave 
functions and for several sets of model parameters || f?J. In the case that the calculated d' 
mass is 600 — 400 MeV above the experimental value 2065 MeV, the decay width T n ~ 16 — 8 
MeV is much too large compared to the experimental result. On the other hand, a decay 
width of rv ~ 10 MeV is still small compared to the width of ~ 150 MeV expected for a 
two-body decay of the d' into the NN* channel with subsequent pionic decay of the N* 
resonance. This channel is not included in the present calculation because it is closed for the 
physically interesting d' masses around 2100 MeV. In the case that the d' mass is close to 
the experimental one, small pionic decay widths are obtained. Using the d' wave functions 
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and parameters of sets 3 & 5 (be's are large), the calculated IVs are both 0.04 MeV which is 
an order of magnitude smaller than the experimental T^ xp . We have also discussed that if one 
could obtain a d' mass of ~ 2100 MeV with a size parameter not much different from that 
of a free-baryon, it would be possible to reproduce the empirical pionic decay width r® xp . 

In this paper the role of the intermediate s 4 p 2 [42] x -configuration in addition to the dom- 
inant s 6 [6]x-one is investigated and found to play an important role in the decay width for 
some parameter sets. The inclusion of the s 4 p 2 [42]x-configuration reduces the magnitude of 
the decay width compared with a calculation using only the s 6 [6]x-configuration. It is also 
shown that the nonstatic correction in the decay operator changes the width even at low- 
momentum transfers and that the form factor at the decay vertex reduces the decay width 
for the parameter sets involving high-momentum transfers. 

We conclude that the experimentally observed small pionic decay width r^ xp ~ 0.5 MeV 
of the d' -dibaryon can be explained by the microscopic shell model wave function and 
the single-quark transition operator provided that the model can reproduce the mass M# 
correctly and the size of d' is not very different from the size of a free baryon. Finally, we 
find it important and necessary to go beyond the plane wave approximation and to solve 
the ttNN final state wave function exactly within the present model in order to obtain an 
improved pionic decay width of the d' -dibaryon. 
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Appendix 

The pionic decay width of the d' -dibaryon is expressed in our model as 

T-jr = 3 r,r- , 



where 
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9max = ^M N (M d , -2M N -m n ) . 



In the above expression the h.o. size parameter b (6/v) denotes that of the d' -dibaryon ( 
free nucleon ). (F.F.) represents the vertex form factor and [NSCL] the nonstatic correction 
effect associated with Eq. (10). 
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